A quantitative study of quasiparticle traps using the single-Cooper-pair-transistor 
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We use radio-frequency reflectometry to measure quasiparticle tunneling rates in the single- 
Cooper-pair-transistor. Devices with and without quasiparticle traps in proximity to the island 
are studied. A 10^ to 10^-fold reduction in the quasiparticle tunneling rate onto the island is 
observed in the case of quasiparticle traps. In the quasiparticle trap samples we also measure a 
commensurate decrease in quasiparticle tunneling rate off the island. 



In a superconductor charge may be transported co- 
herently by Cooper-pairs and incoherently by quasiparti- 
cles. For nanoscale Coulomb blockade devices such as the 
single-Cooper-pair-transistor (SCPT) and Cooper-pair- 
box (CPB), these two modes of transport are mutually 
exclusive 0, 0, 0] • Coherent charge transport is tem- 
porarily halted while a quasiparticle is present on the 
device island. This effect, caused by the electrostatic 
energy cost of the quasiparticle, is known as quasiparti- 
cle poisoning and leads to Cooper-pair tunneling being 
stochastically switched on and off. Quasiparticle poison- 
ing is usually unwanted, especially in charge-based super- 
conducting qubits (e.g. the CPB) where it changes the 
qubit bias point in a random telegraphic way and places 
limits on quantum coherence [1, 0, Q • 

In order to suppress quasiparticle poisoning the density 
of quasiparticles in the device leads should be reduced. 
In principle this can be simply achieved by lowering the 
temperature, however this is not always sufficient and 
an additional approach is to use so-called "quasiparticle 
traps". These are normal metal regions contacting the 
superconducting leads in proximity to the device island 
[H. In this paper we measure quasiparticle tunneling 
rates for SCPTs with and without quasiparticle traps. 
We report quantitatively on the effect of the traps, and 
discuss the implications of our results for designing su- 
perconducting Coulomb blockade devices with reduced 
quasiparticle poisoning. 

Quasiparticle traps provide a sink to reduce quasipar- 
ticle density. When a superconductor is brought in con- 
tact with a normal metal, an induced superconducting 
gap develops in the normal metal, this is known as the 
proximity effect. Quasiparticles in the superconductor 
diffuse into the normal metal side where they relax to- 
wards the gap edge on a timescale which depends on the 
electron-phonon coupling The trapped quasiparti- 

cles eventually recombine, emitting a phonon with insuf- 
ficient energy to break a Cooper pair in the superconduc- 
tor. Quasiparticle traps have been studied in the context 
of x-ray photon detection and cooling through 



poisoning in SCPTs [H, copper leads were contacted at a 
distance of <1 jam from the junctions and a 2e-periodic 
super current was observed at low temperatures. Quasi- 
particle traps have often been used where quantum co- 
herence is of interest |5|, [III, Q, [H, [l^- A variety of 
results have been achieved with immunity from quasi- 
particle poisoning not always evident [l^, 17 1: therefore 
providing a motivation for this work. 

We make use of the enhanced superconducting gap of 
aluminum with decreasing film thickness [18]. The SCPT 
island gap is made greater than the lead gap (A^ > A^) 
to reduce the depth of the quasiparticle potential well on 
the island [2]. Also, the leads of the SCPTs are isolated 
from the normal metal contacts by a region of higher su- 
perconducting gap, forming a quasiparticle barrier. This 
was designed to reduce the effect of quasiparticles in the 
leads being trapped by the normal metal bond-pads, al- 
lowing the full effect of the intentional traps to be ob- 
served [Fig. [T](c)]. The island and quasiparticle barrier 
layers are 7 nm thick (A^ = 298 ± 9.4 jaeV) while the 
leads are 30 nm thick (A/ = 209 ± 11 /ieV) De- 
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first applied quasiparticle traps to reduce quasiparticle 



FIG. 1: (a) Micrograph of a QT device with an island volume 
voli—7^0 nmxl25 nmx7 nm. Bright features are the nor- 
mal metal traps 200-300 nm from SCPT junctions, (b) The 
rf- circuit has a resonance frequency of ^ 320 MHz and in- 
cludes a chip inductor (L = 470 nH), a parasitic capacitance 
{Cp = 0.53 pF) and the SGPT. (c) Energy gap profile showing 
isolation of the leads from the normal metal bond-pads. 
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FIG. 2: (a) Histogram of times spent in the even and odd 
states for a QT device, solid lines are fitted exponentials. Here 
todd = 3.1 ms and teven — 964 jjLS. Inset: switching between 
even and odd states at = 1. (b) Upper: Averaged super- 
current oscillations and Lower: energy band diagram of the 
SCPT at zero current bias. The quasiparticle bands (dashed 
Une) are lifted by an energy creating a potential well of 
energy 5E relative to the lower Cooper pair band (solid line), 
(c) Allowed transitions of a quasiparticle including the relax- 
ation rate to the bottom of the potential well on the island. 



vices were fabricated with and without direct contact to 
quasiparticle traps. We refer to these device types as 
QT (quasiparticle traps) and NT (no traps) respectively 
[Fig. [TJa)]. AuPd alloy was chosen as the trap metal 
as we find this reliably gives a low contact resistance to 
aluminum. Our quasiparticle traps were typically 200- 
300 nm from the SCPT tunnel junctions and have ap- 
proximately the same volume as the leads {vol = 2x10 
/imxlOO nmx30 nm). 

The SCPT bandstructure determines the energetics of 
quasiparticle tunneling [Fig. [2](b)]. The difference in su- 
perconducting gaps between the leads and the island lifts 
the energy minima of the quasiparticle state by an en- 
ergy = A^ — A/ In our devices, the quasiparticle 
state remains lower in energy than the ground band of 
the SCPT at = 1. Therefore a quasiparticle potential 
well is formed on the SCPT island with a depth given by 
5E = Ec— ^ —SA [2] . The presence of a well means that 
a quasiparticle in the leads can tunnel into the island, re- 
lax, and after some time can be thermally activated out. 
We refer to the time constants for the poisoning process 
as teven (the cvcu state lifetime) and the unpoisoning pro- 
cess as todd (the odd state lifetime). From normal state 
Coulomb diamonds {B = 3T) we find the charging en- 
ergy of the principal QT (NT) device as, Ec = e'^/2C^ = 
170 fieV (185 /ieV); the total device resistance at 4.2 K 
is 42 (54 kft); and the Josephson energy per junction 



from the Ambegoakar-Baratoff relation is Ej = 37 /ieV 
(28 /ieV). The SA we expect from SIS junction measure- 
ments is 89 ± 12 /ieV. Hence, the expected depths of 
the island potential well are SE = 62.5 jneY (82 //eV), in 
both cases SE ^ kT. For our set of devices Ec and Ej 
are consistent to within 25%. 

Radio-frequency reflectometry is employed to perform 
high-bandwidth measurements of quasiparticle tunneling 
[l9|. This consists of embedding the SCPT in a reso- 
nant circuit and detecting the amplitude (and phase) of 
a small rf signal reflected from the circuit at resonance 
[Fig. [TJb)]. The reflected carrier signal passes through a 
low noise cryogenic amplifier before being further ampli- 
fied and demodulated. Briefly, the presence of a quasi- 
particle on the device island shifts the supercurrent oscil- 
lation by a gate charge of e. This changes the impedance 
of the SCPT and changes the amplitude and phase of the 
reflection coefficient [3[Q. 

We monitor the demodulated signal rf signal (Vqut) 
with the device biased on a supercurrent peak (upper 
graph Fig. El^b)) and observe two- level switching due to 
quasiparticle tunneling [Fig. EJa) inset]. Traces consist 
of 10^ data points and ~ 1 x 10^ switching events, and 
are analyzed by comparing Vqut to a median value. We 
plot a histogram of times spent in each state [Fig. [2ja)], 
an exponential is fitted to the histogram and the de- 
cay constant defines the state lifetimes. A good fit in- 
dicates Poissonian tunneling processes. Finite receiver 
bandwidth causes a systematic overestimate of time con- 
stants, and we follow to correct for this. Low pass 
filters of between 1.0 and 10.7 MHz are used and the 
receiver bandwidth is estimated for each filter. We op- 
erate at an rf-bias where the time constants are not sig- 
nificantly perturbed, corresponding to a carrier power of 
-100.5 dBm (-104 dBm) for the QT (NT) devices [3]. 

For the QT device a dramatic increase in tg^en (re- 
duction in poisoning rate), compared to the NT de- 
vices, is observed across the whole temperature range 
[Fig. [31(a)]. The saturated low temperature values are 
teven = 1 uis for the QT case and teven = 1-9 /is for the 
NT case. The equilibrium expression for quasiparticle 
density is riqp = D{eF)V^7rAikT exp(— A///cT), where 
D{eF) = 1.45x10^'^ m~^J~^ is the normal-state den- 
sity of states (including spin) of aluminum at the Fermi 
energy. Since t~J^^^ is proportional to the quasiparticle 
density in the leads, we fit to the following expression 
^even = ^ {^qpi^qp, A/) + nqp{T, A/)), whcrc B IS 8i Con- 
stant, Tqp is the saturation temperature and Uqp refers to 
the quasiparticle density in the leads \3] . For the NT and 
QT devices, we find the minimum quasiparticle tempera- 
tures are Tqp = 248 ± 3 mK and Tqp = 204 ± 1 mK; the 
superconducting gaps Ai = 210 ± 20 /ieV and A^ = 120 
± 3 /ieV; and the constant B = 2.2 ± 1.9x10"^^ m^s"^ 
and B = 4.0 ± 0.6x10"^^ m^s-^. 

A reduction in quasiparticle temperature is seen for 
the QT device, however this only partially explains the 
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FIG. 3: Temperature dependence of (a) teven and (b) todd for 
a QT device (triangles) and a NT device (squares). Solid lines 
are fits to the quasiparticle density in the leads (see text). 



change in tg^en- If the effect of the quasiparticle traps 
was purely to thermalize the leads we would expect the 
temperature dependencies to overlie one another at high 
temperature with the QT case saturating at a lower tem- 
perature. To illustrate this further we take the expression 
for UqpiT) and estimate the difference in teven for the 
observed change in minimum quasiparticle temperature, 
n,p(248 mK)/n,p(204 mK) = te^en(204 mK)/te^en(248 
mK) ~ 9: a relatively small fraction of the total change. 

The NT gap parameter extracted for the leads (A/ = 
210 ± 20 /ieV) agrees with our expected value, but a sig- 
nificant suppression is seen for the QT device (A^ = 120 
± 3 /ieV). We check if this corresponds to a real reduc- 
tion in A/ by performing a voltage-biased measurement 
of quasiparticle tunneling threshold. This threshold oc- 
curs at Vds =2eA^+2eAi = 998 /ieV for the QT device 
so, assuming an unchanged eA^ =298 jieV then eA^ = 
201 jieV which is close to the NT case. We conclude that 
that A^ is not significantly changed by the presence of the 
normal metal. However, it is not the only quasiparticle 
activation energy in a superconductor- normal bilayer. 

Furthermore we can look at the values of the propor- 
tionality constant, 5, for the tunneling rate. Theoret- 
ical analysis of the tunneling rates show that t~^en — 



Gn 



[2l m, where Vi,i{Ek) = 



is the quasiparticle density of states and G 



Gi + ^2, in S, is the total tunnel barrier conductance. 
Thus our proportionality constant is predicted to be 5 = 
1.6x 10~^^ m^s~^ showing close agreement with the mea- 
sured value in the NT case but not in the QT case. 
From these fitting parameters it appears that the den- 
sity of quasiparticle excitations, and their rate of tun- 
neling onto the SCPT island, in the NT case is well- 
described. However, from the QT device, we observe that 
the quasiparticle density in the superconducting side of 



a superconductor-normal bilayer is significantly reduced 
and is not given by the equilibrium expression for quasi- 
particle density nqp{Teff) at an effective temperature 
Teff. We note that the increased teven for the QT devices 
is reproduced amongst our sample set which consists of 
3 of each device type [Fig. [31(a)]. 

We now discuss qualitatively how the quasiparticle 
density is affected by the presence of a quasiparticle trap. 
In a superconductor the number of quasiparticles {Nqp) 
and phonons (N^) with E > 2 A are closely c oup led and 
described by the Taylor- Rothwarf equations [231] . In the 
presence of a quasiparticle trap, it is necessary to add ad- 
ditional terms to describe loss to {-TrNqp) and gain from 
(TuN^p) the trap [23]. Quasiparticle population of the 
trap is given by Nj^. The constants G and R represent 
the quasiparticle generation and recombination rates in 
the superconductor. The quasiparticle master equation 
for the superconducting film follows. Other equations 
treat the quasiparticle density in the trap (A^^) and the 
phonon density in both films. 



qp 



dt 



2GN^ 



vol 



(1) 



To fully explain the teven behavior of the QT sam- 
ples, the solution of this equation should be fitted to the 
data in Fig. [31(a). We do not perform this analysis but 
note that the significant reduction in quasiparticle den- 
sity of the superconducting film implies that trapping, 
rather than recombination becomes the dominant loss 
term. In this case, where tunneling to the trap is the 
fastest timescale, the quasiparticle density in the lead is 
given by detailed balance between the lead and the trap, 
^rNqp = TijNqp. The trapping (F^) and un-trapping 
rates (Tu) will depend on the lead-trap interface. Since 
Nqp is smaller for QT devices, it follows that the number 
of recombination phonons (A^a;) is also reduced. 

Now we turn to the behavior of todd- Two of three 
QT samples show a significant increase in todd over the 
NT devices [25]. The increase is comparable to that for 
teven- there is a factor of 9000 for the principal QT and 
NT devices [Fig. [31(b)]. We can estimate the reduction in 
temperature needed from oc exp {—dE/kT). Taking 
the mean value of quasiparticle potential depth {dE = 72 
fieV) for the two principal devices, a reduction in phonon 
temperature from 248 mK to 67 mK can explain the dif- 
ference observed in todd- Phonons (with E > 2A), re- 
sulting from quasiparticle recombination should be con- 
sidered in addition to the thermal phonons. If the effect 
of the thermal phonons was to dominate the un-poisoning 
process, then todd should saturate at the phonon temper- 
ature of the sample. We expect the phonon tempera- 
ture to be lower than the quasiparticle temperature and 
close to the base temperature of the dilution refrigera- 
tor (T ~ 50 mK). In these devices todd saturates at a 
temperature close to teven implying that recombination 
phonons may play a significant role in the unpoisoning 
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FIG. 4: (a) Values of todd and teven plotted for QT (trian- 
gles) and NT devices (squares) at fridge base temperature, 
(b) The probability of having a quasiparticle on the island 
(Podd), determined from the time constants at this tempera- 
ture. QT(NT) devices are indicated as filled (unfilled). The 
principal QT(NT) device is noted with *. 



process. The relative effect of thermal and recombina- 
tion phonons may depend on sample parameters. In a 
previous device, with lower charging energy [3], todd was 
observed to saturate at a lower temperature that teven- 

In the QT device with the short todd = 880 ns 
[Fig. SJa)] we see no temperature dependence of todd^ 
suggesting the quasiparticle exits without thermalization. 
This behavior may be explained by the existence of two 
unpoisoning mechanisms [21]. Either the quasiparticle 
tunnels off the island immediately in an elastic unpoi- 
soning process [Fig.[2jc)]. Or, it relaxes to the bottom of 
the potential well with a rate 7 which strongly depends 
on the island potential depth (7 oc dE^-^). Subsequently 
thermal activation is required for the unpoisoning event 



to occur. The escape rate todd 



5E 



for 



2e2 ^ V-^ ) 

the first scenario is governed by the conductance of the 
junctions, the quasiparticle level spacing on the island, 
and the initial energy of the quasiparticle (of order 
T) [2l|, We calculate todd = 280 ns, in reasonable 
agreement with the measured value. The reason for the 
dominance of the elastic process in this sample is not 
clear, one possibility is variation of the thin-films leading 
to a reduced potential well depth on the island. 

Our measurements are sensitive to unpoisoning and 
poisoning processes, however this is not always the case. 
For example, in electrometry of a CPB by a single- 
electron-transistor, averaged measurements of the box 
charge are typically performed 0. This yields the av- 
erage probability of having a quasiparticle on the island 
Podd without directly providing details of the individual 
tunneling rates. Taking our measured times for the QT 
(NT) device then podd = , — = 0.75 (0.15) at base 
temperature. The time the SCPT remains poisoned is 
greater in the QT device, even though the the tunneling 
rate onto the island is significantly reduced. Thus one 
might conclude from Podd that the quasiparticle traps 
actually had a negative effect. Figure 4(b) plots Podd for 
all the samples measured and shows no obvious trend be- 
tween the QT and NT devices. This is consistent with 



the electrometry measurements of a CPB observing no 
clear benefit of quasiparticle traps fid\. 

In conclusion, quasiparticle traps can definitely help re- 
duce quasiparticle poisoning in superconducting qubits: 
a reduction of 2-3 orders of magnitude in poisoning rate 
was observed here. However, quasiparticle traps are 
only part of the solution to the quasiparticle poison- 
ing problem. They should be used in conjunction with 
small charging energy devices and graded superconduct- 
ing gaps to ensure high unpoisoning rates. 
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